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We present measurements of B"* 



D*" 



and B^ 



U 7 



decays in a data sample of 



657 X lO'' BB pairs collected with the Belle detector at the KEKB asymmetric-energy e"'"e collider. 
We find 446^56 events of the decay B^ — )• D'^W^Vr with a significance of 8.1 standard deviations, 
and 146^41 events of the decay B^ — >■ D'^r'^i/r with a significance of 3.5 standard deviations. The 
latter signal provides the first evidence for this decay mode. The measured branching fractions are 
B{B+ -^ D'°T+Vr) = (2.12t[!-^*(stat) ±0.29(syst))% and B(B+ -> D°T+Ur) = (0.77 ± 0.22(stat) ± 
0.12(syst))%. 



PACS numbers: 13. 20. He, 14.40. Nd 



Measurements of leptonic and semileptonic decays of 
B mesons to the r lepton can provide important con- 
straints on the Standard Model (SM) and its extensions. 
Due to the large mass of the lepton in the final state, 
these decays are sensitive probes of models with extended 
Higgs sectors [l|]. Semileptonic modes with b — )■ CT~Dr 
[2] transitions provide more observables sensitive to new 
physics than purely leptonic B'^ — > t~^Vt decays. Of 
particular interest is r polarization. The effects of new 
physics are expected to be larger in B — > Dt^v^- than in 
B — >■ D*t'^Vt. We note that decays to the vector meson 
offer the interesting possibility of studying correlations 
between the D* polarization and other observables [3[. 

The predicted branching fractions, based on the SM, 
arc around 1.4% and 0.7% for 5° -J> D*~T^Vr and 
B^ — > D^T^v^, respectively (see e.g., [J|). Despite rela- 
tively large branching fractions, multiple neutrinos in the 
final states make the search for semi-tauonic B decays 
very challenging. Inclusive and semi-inclusive branch- 
ing fractions have been measured in LEP experiments [5| 
with an average branching fraction of B{h — )■ tvtX) = 
(2.48±0.26)% [6|. The exclusive decay was first observed 
by Belle [7|] in the B° — > D*~T^Vr mode. Other modes 
have also been measured by BaBar [8| and Belle [9| . The 
results are still statistically limited. In particular, the 
Belle prelimnary result [9| is the only evidence to date 
for B'^ — > D'^t'^Vt-. Further improvements in precision 
could tightly constrain theoretical models. 

Decays of B mesons to multi-neutrino final states can 
be studied at i3-factories via the recoil of the accom- 



panying B meson (Stag)- Reconstruction of the Stag 
allows one to calculate the missing four-momentum in 
the Bsig decay; this helps separate signal events from co- 
pious backgrounds. At the same time the presence of 
a Stag strongly suppresses the combinatorial and con- 
tinuum backgrounds. The disadvantage is the low Stag 
reconstruction efficiency. To increase statistics, we re- 
construct the Stag "inclusively" from all the remaining 
particles after the Sgjg selection (see Ref. [7|). A data 
sample consisting of 657 x 10® BB pairs is used in this 
analysis. It was collected with the Belle detector [10| at 
the KEKB asymmetric-energy e+e" (3.5 on 8 GeV) col- 
lider [ll| operating at the T(4S') resonance (-y/s = 10.58 
GeV). 

We use Monte Carlo (MC) simulations to estimate 
signal efficiencies and background contributions. Large 
samples of the signal S+ —J- D^*'^^t^Vt decays are gen- 
erated with the EvtGen package [l^ using the ISGW2 
model [l3[ . Radiative effects are modeled using the PHO- 
TOS code [ij] . We use large MC samples of continuum 



qq [q = u, d, s, c) and inclusive BB events to model the 
background. The sizes of these samples are, respectively, 
six and nine times that of the data. 

Primary charged tracks are required to have impact 
parameters consistent with an origin at the interaction 
point (IP), and to have momenta above 50 MeV/c in 
the laboratory frame. Kg mesons are reconstructed us- 
ing pairs of charged tracks satisfying 482 MeV/c^ < 
M^+^- < 514 MeV/c^ with a vertex displacement from 
the IP consistent with the reconstructed momentum vec- 



tor. Muons, electrons, charged pions, kaons and protons 
are identified us ing information from particle identifica- 
tion subsystems [15[ . The momenta of particles identified 
as electrons are corrected for bremsstrahlung by adding 
photons within a 50 mrad cone along the lepton trajec- 
tory. 

The TT*^ candidates are reconstructed from photon pairs 
having 118 MeV/c^ < M^-y < 150 MeV/c^. For candi- 
dates that share a common 7, we select the one with the 
smallest x^ value resulting from a tt^ mass-constrained 
fit. To reduce the combinatorial background, we require 
that the photons from the 7r° have energies greater than 
60 McV - 120 MeV, depending on the photon's polar 
angle. Photons that are not associated with a tt" are 
accepted if their energy exceeds a polar-angle dependent 
threshold ranging from 100 MeV to 200 MeV. 

The Z)° candidates are reconstructed in the A'+tt^ and 
A'+tt^tt'^ final states. We accept D^ candidates having 
an invariant mass in a 3cr window of the nominal il/jjo 
mass. 

The D*^ candidates are reconstructed from D^ir'^ . We 
require that the mass difference AM = Md* —Mj^o is in a 
3(7 window around its nominal value. We also accept D^^ 
pairs that do not fulfill the requirement on AM if they are 
kinematically consistent with the hypothesis that Z)° and 
7 come from the decay Z)*° — ?► Z)°7r° with one undetected 
photon (7niiss) from the 7r° decay ("partial reconstruc- 
tion" of Z)*°). For this purpose cos(0^,-y^;^J, the cosine 
of the angle between two photons from the tt" is calcu- 
lated in the D*^ rest frame taking the nominal D*^ and 
7r° masses. We require | cos(6'^,-y,^j^J| < 1.1 (taking into 
account experimental precision) and that the energy of 
the detected photon exceeds 120 MeV. The partial recon- 
struction of Z)*° increases the reconstruction efficiency by 
a factor of about four, but due to higher background it 
is only used in the subchannels with D° — > A'+tt^ decay. 

To reconstruct the r lepton candidates, we use the 



-> /i+i'^ inl- 



and 



modes. 



In the latter case, we also take into account the contri- 
bution from the r"'' — > p^Dr channel. The t+ — > tt^Dt 
mode has a sensitivity similar to the t+ — > e^v^DT- or 
T+ — > n^v^Vr mode, and can be used to study t polar- 
ization. For this channel, due to the higher combinatorial 
background, we analyze only the decay chains with the 
D^ — > A'+TT^ mode. In total, we consider 13 different de- 
cay chains, eight with D*^ and five with Z)" in the final 
states. 

The signal candidates are selected by combining a 
]j{*)0 nieson with an appropriately charged electron, 
muon or pion. In the sub-channels with the r"*" — >■ tt^Dt 
decay, the large combinatorial background is suppressed 
by requiring the pion energy E^^ > 0.6 GeV. From mul- 
tiple candidates we select a (i)^*)°ci+) pair (throughout 
the paper dr stands for the charged r daughter: e, /i or 
tt) with the best i)^*^° candidate, based on the value of 
AM (for subchannels where AM is available) or Mjjo. 



For the pairs sharing the same !)(*)'' candidate, we select 
the candidate with the largest vertex probability fit on 
the tagging side. 

Once a Bgig candidate is found, the remaining par- 
ticles that are not assigned to _Bsig are used to recon- 
struct the -Btag decay. The consistency of a -Stag candi- 
date with a i?-meson decay is checked using the beam- 
energy constrained mass and the energy difference vari- 
ables: Mtag = y^£^bcam " P?ag> Ptag = T,^P^, and 
AEtag = Etag — Abeam, ^tag = J2i ^i^ whcrC Abeam IS 

the beam energy and Pi and Ei denote the 3-momcntum 
vector and energy of the i'th particle. All quantities are 
evaluated in the T(4S') rest frame. The summation is 
over all particles that are assigned to Stag- We require 
that the candidate events have Mtag > 5.2 GeV/c^ and 
-0.3 GeV < AAtag < 0.05 GcV. With this require- 
ment the Mtag distribution of the signal peaks at the B^ 
mass with about 80% of the events being contained in 
the signal-enhanced region Aftag > 5.26 GeV/c^. 

To suppress background and improve the quality of 
the Btag selection, we impose the following requirements: 
zero total event charge; no charged leptons in the event 
(except those coming from the signal side); zero net pro- 
ton/ antiproton number; residual energy in the electro- 
magnetic calorimeter (i.e., the sum of energies that are 
not included in the i?sig nor Stag) should be less than 
0.35 GcV (0.30 GeV or 0.25 GeV in sub-channels with 
higher backgrounds); the number of neutral particles on 
the tagging side iV^o -I- A^^ < 6, A^^ < 3, and less than 
four tracks that do not satisfy the requirements imposed 
on the impact parameters. For decay modes with higher 
background, we impose further constraints on the total 
event strangeness and require no K^ in the event. These 
criteria, which we refer to as "the Stag-selection", re- 
ject events in which some particles were undetected and 
suppress events with a large number of spurious show- 
ers. In the samples of the (Z)(*^°Z"'") pairs {I = e,fj,), the 
dominant background comes from semileptonic B decays, 
B+ -^ D'^*'>°Xl+iyi, whereas in the case of the (Z)(*)°7r+) 
pairs, the combinatorial background from hadronic B de- 
cays dominates. 

Further background suppression exploits observables 
that characterize the signal decay: missing energy 
Emiss = Abeam-£^_D(.)o-A^+; visiblc energy Avis , «-e., the 
sum of the energies of all particles in the event; the square 



of missing mass M, 



2 

miss 



E' 



(Psig - PB(*)o - P^+) 



and the effective mass of the (r+i^T-) pair, q^ = (Abe 
EeH')°Y - (Psig - Pnoo)^ where Psig = -Ptag (all kine- 
matical variables are in the T(4S') rest frame). The most 
useful variable for separating signal and background is 
obtained by combining Amiss and (Z)(*)°d+) pair momen- 



tum: ATmiss = (Amiss - |Pj5(*)0 + Pd+D/y^bcam " "^5 + 

where mg+ is the A"*" mass. The A"iniss variable is closely 
related to the missing mass in the Asig decay but does 



not depend on the -Stag reconstruction [7[. 

The signal selection criteria are optimized individu- 
ally in each decay chain, by maximizing the expected 
significance Ns/y/^S + ^b, where Ns and Nb are the 
number of signal and background events in the signal- 
enhanced region, assuming the SM prediction [j] for the 
signal branching fractions. The expected background 
Nb is evaluated using generic MC samples. We require 
Svis < 8.3 GeV - 8.5 GeV, E^^iss > 1-5 GeV - 1.9 GeV 
and Xmiss > 2.0 - 2.75 for leptonic r decays, or Xmiss > 
1.0 - 1.5 for the modes with r — >■ ttj/t-. In the latter case, 
where the r decays to a final state with a single neutrino, 
we further require cos{9,y-^,j^) to be in the range [—1,1], 
where 9u^iJ2 denotes the angle between the two neutri- 
nos in the [t^Vt) rest frame and is calculated from the 
M^jgg and q^ variables. In the sample with [D^d^) pairs, 
to suppress the cross-feeds from the B —J- D*T'^iyr modes, 
we impose a loose requirement on q^ < 9.5 GcV^/c"'. 

The above requirements result in flat Mtag distribu- 
tions for most background components, while the distri- 
bution of the signal modes remains unchanged. The main 
sources of the peaking background are the semileptonic 
decays B+ -> D*°l+vi and B -)■ D'-*\l+vi (including 
D**l+vi). 

In order to estimate the peaking background reliably, 
in particular from poorly known semileptonic modes of 
the type B -^ D**lvi, wc divide the MC sample into 
the following categories: B — )■ D*l'^i/i, B -^ Dl~^i/i, 
B — >■ D**l~^i'i, other B decays, cc and {uu + dd + ss) 
continuum. The normalizations of these components are 
determined from simultaneous fits to experimental dis- 
tributions of Mtag, A£;tag, Ed^, Xmiss, £^vis, 9^, and 
i?2 , the ratio of the second and zeroth Fox- Wolfram mo- 
ments [16| . These fits are performed separately for the 
subsamples defined by the (I)'*^°(i+) pairs, excluding the 
region Mtag > 5.26 GcV/c and Xmiss > 2.0, where wc 
expect enhanced signal contributions. 

The signal and combinatorial background yields are ex- 
tracted from an extended unbinned maximum likelihood 
fit to the Mtag and Pjjo (momentum of D^ from B^ig mea- 
sured in the T{AS) frame) variables. The Aftag variable 
allows us to separate the combinatorial background from 
the signal, while Pqo helps to distinguish between the 
two signal modes. Correlations between these variables 
are found to be small. 

Paramctcrizations of two-dimensional probability den- 
sity functions (PDFs) are determined from the MC 
samples. They are expressed as the product of onc- 
dimcnsional PDF's for each variable. The PDF's for Mtag 
of the signal and peaking background components are de- 
scribed using an empirical parameterization introduced 
by the Crystal Ball collaboration [17|, while combinato- 
rial backgrounds are parameterized by the ARGUS func- 
tion [18|. It has been empirically found that the PDF's 
for P]jo are well modeled as a sum of two Gaussian dis- 
tributions. 



The fits are performed in the range Mtag > 5.2 
GeV/c^, simultaneously to all data subsets. In each of 
the subchannels, we describe the data as the sum of four 
components: signal, cross-feed between D*^t~^i't and 
D^t~^Vt, combinatorial and peaking backgrounds. The 
common signal branching fractions B{B~^ — >■ D*^T^Vr) 
and B{B^ — > D^t'^Vt), and the numbers of combinato- 
rial background in each subchannel are free parameters of 
the fit, while the normalizations of peaking background 
contributions are fixed to the values obtained from the 
rescaled MC samples (as described above). The signal 
yields and branching fractions for _B+ -^ D''*^^t^i't de- 
cays are related using the following formula, which as- 
sumes equal fractions of charged and neutral B meson 
pairs produced in T(45) decays: B{B+ -^ D'^*'^^T+Vr) ^ 
N^ * /{Nbb X J2k ^kBk), where Nbb is the number of 
BB pairs and the index k runs over the 13 decay chains; 
ek denotes the reconstruction efficiency of the specific 
subchannel and Bk is the product of intermediate branch- 
ing fractions. All the intermediate branching fractions 
arc taken from the PDG compilation [6[. The efficiencies 
of the signal reconstruction, as well as the expected com- 
binatorial and peaking backgrounds are given in Table ID 

The signal extraction procedure has been tested by 
fitting ensembles of simulated experiments containing 
all signal and background components. These pseudo- 
experiments are generated using the shapes of the fit- 
ted PDF's for the signal and background components 
and with the number of events are Poisson-distributed 
around the expected yields. The pull distributions of 
the extracted signal branching fractions are consistent 
with standard normal distributions. The small biases 
in the mean values are included in the final systematic 
uncertainties. The procedure established above is ap- 
plied to the data. The Mtag and Pdo distributions for 
the Z)*°r+i^r and D'^r^Vr samples in data are shown in 
Fig. [TJ The overlaid histograms represent the expected 
background, scaled to the data luminosity. A clear ex- 
cess of events over background is visible in the signal- 
enhanced region. 

The branching fractions extracted from the fit are 
B{B+ -^ D*^T+Vr) -= (2.12+°-28(stat))% and i3(B+ ^ 
D°T+Vr) = (0.77 ± 0.22(stat))%. The signal yields 
are 4461^^ B 



D*°T+Vr events and U6tfi B+ 



D^T^Vr events. The statistical significances, defined as 
S — ^— 21n(£o/>Cmax), corespond to 8.8 and 3.6 stan- 
dard deviations {a), respectively. Here £max denotes the 
maximum likelihood value and Cq is the likelihood for 
the zero signal hypothesis. The fitted yields of combina- 
torial background in the individual submodes are consis- 
tent within statistical uncertainties with the MC-based 
expectations. The fit results are summarized in Table HI 
The fit projections in Mtag and Ppo are shown in Fig. [T] 

As a cross-check, we extract the signal yields from 
an extended unbinned maximum likelihood fit to one- 



TABLE I. The yields of signal (Ns) and combinatorial background (Nb) events determined from fits to data, number of expected 
combinatorial (N^) and peaking (Np) background events, signal selection efficiencies (e = X^feefe-Sfe), extracted branching 
fractions (B) and statistical significances (E). A'j, and A^^ are evaluated from fits of the generic MC samples to experimental 
distributions. The numbers in parentheses refer to the signal D*^{D") modes reconstructed as D°{D*^^). The efficiencies include 
intermediate branching fractions. The listed errors are statistical only. The results are summed over the considered _D'*'° and 
T decay modes. 



Mode 



Ns 



N„ 



N, 



JMT 



N„ 



e(10- 



B{%) E(a) 



D*°T+u^ 446+56(226) 1075+^' 
i46+*2 



, + 0.1^8 



DV-t 



?(15) 



1245t« 



1029 ± 20 31.0 ± 17.7 32.6 ± 0.2(16.3) 2.121o 27 
1310 ± 19 78.2 ± 12.6 30.0 ± 0.4(3.2) 0.77 ± 0.22 3.6 
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FIG. 1. The fit projections to Mtag, and Pqo for Mtag > 
5.26 GeV/c^ (a,b) for D*°T+Vr, (c,d) for D°T+iyr- The black 
curves show the result of the fits. The solid dashed curves 
represent the background and the dashed dotted ones show 
the combinatorial component. The dot-long-dashed and dot- 
short-dashed curvess represent, respectively, the signal con- 
tributions from B+ -> D*°T+Ur and B+ -^ D°T+Vr- The 
histograms represent the MC-predicted background. 



of real processes. The uncertainty in the i?tag and part of 
the _Bsig reconstruction efficiency is evaluated from data 
control samples with B+ -> i)*°7r+ and 5+ -^ i)°7r+ de- 
cays on the signal-side. The absolute normalizations of 
the data and MC control samples agree to within 13%. 
The difference, as well as uncertainties in the relative 
amounts of D*" — D° cross-feeds are included in the sys- 
tematic uncertainty of -Btag and B^ig reconstruction. The 
remaining uncertainties in the lepton identification and 
signal selection are estimated separately. The latter are 
determined by comparing MC and data distributions in 
the variables used for signal selection. The uncertainties 
due to the partial branching fractions Bk are taken from 
the errors quoted by the PDG [6|. 

The systematic uncertainties in the signal yield origi- 
nate from the background evaluation and from the PDF 
parameterizations of the signal and background compo- 
nents. The resulting error is evaluated from changes in 
the signal yields obtained from fits where the PDF pa- 
rameters and the relative contributions of the background 
components are varied by zblcr. 

All of the above sources of systematic uncertainties 
arc combined together taking into account correlations 
between different decay chains. The combined systematic 
uncertainty is 13.9% for the B'^ — >■ D*'^t'^Vt mode and 
15.2% for B+ -^ D°T+Vr. 



dimensional distributions in Mtag and obtain consistent 
results with the two-dimensional fit. We also examine 
the distributions of variables used in the signal selection, 
applying all requirements except those that are related 
to the considered variable. In all cases the distributions 
are well-reproduced by the sum of signal and background 
components with normalizations fixed from the fit to the 
(Mtag, -Pdo) distribution. 

The systematic uncertainties in the branching fractions 
are summarized in Table |lll They include uncertainties 
in the total number of BB pairs, the effective efficiencies 
J2k^kBk^ and the signal-yield extractions. The system- 
atic uncertainties associated with the effective efficien- 
cies include errors in determination of the efficiencies for 
Stag reconstruction and (-D(*)°(i+) pair selection, coming 
from efficiencies of tracking, neutral particle reconstruc- 
tion, particle identification, and from imperfect modeling 



TABLE II. Summary of the systematic uncertainties. 



Source 


D*V+^. 


D'^r+u^ 


Nbb 


±1.4% 


±1.4% 


Reconstruction of Btag and Bsig 


±12.9% 


±12.8% 


Lepton-id and signal selection 


+ 1.5(v 
-1.6 /C 


+4.4q/ 


Shape of the signal PDF's 


±2.5% 


±6.0% 


Comb, and peaking backgrounds 


±3.3% 


±2.7% 


Fitting procedure 


±0.8% 


±1.5% 


Total 


±13.9% 


±15.2% 



We include the effect of systematic uncertainties in the 
signal yields on the significances of the observed signals 
by convolving the likelihood function from the fit with a 
Gaussian systematic error distribution. The significances 
of the observed signals after including systematic uncer- 
tainties are S.lcr and 3.5cr for the 5+ -^ D*^t^Vt and 
B^ — ^ D^t^Vt modes, respectively 



In conclusion, in a sample of 657x10^ 
measure branching fractions B{B^ -^ 



(2.12t°J8(stat)±0.29(syst))%, and S(S^ 



BB pairs wc 

(0.77 ± 0.22(stat) ± 0.12(syst))%, which are consistent 
within experimental uncertainties with SM expectations 
[J|. The result on i?+ — > D^r^v-r is the first evidence for 
this decay mode. 
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